The zero age main sequence of WIMP burners 
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We modify a stellar structure code to estimate the effect upon the main sequence of the accretion 
of weakly interacting dark matter onto stars and its subsequent annihilation. The effect upon the 
stars depends upon whether the energy generation rate from dark matter annihilation is large enough 
to shut off the nuclear burning in the star. Main sequence WIMP burners look much like protostars 
moving on the Hayashi track, although they are in principle completely stable. We make some brief 
comments about where such stars could be found, how they might be observed and more detailed 
simulations which are currently in progress. Finally we comment on whether or not it is possible to 
link the paradoxically young OB stars found at the galactic centre with WIMP burners. 



There is growing evidence that for each gram of bary- 
onic matter in the universe, there are around five grams 
of dark matter which does not couple to the electromag- 
netic force fl| . One of the more convincing candidates for 
this dark matter are weakly interacting massive particles 
(WIMPs). Because WIMPs have masses and couplings 
close to the weak scale, they naturally give rise to a relic 
abundance of dark matter close to that observed after 
freeze out in the early universe. 

Because WIMPs couple weakly to standard model 
particles, there is a small but non-zero WIMP-nucleon 
cross section. The tightest constraints for the spin- 
independent cross section come from the XENON ex- 
periment ■ The WIMP-nucleon interaction means that 
some WIMPs are gravitationally captured by stars, a pro- 
cess which has been well studied 3, 4] . Usually the focus 
of such investigations is the possibility of dark matter an- 
nihilating into high energy neutrinos, which then escape 
the star to be potentially detected by neutrino experi- 
ments like IceCube. If the accretion of dark matter were 
large enough however, one might expect that the stars 
themselves could change. 

The first way WIMPs can affect a star is by annihilat- 
ing into standard model particles in its core, providing 
another energy source in addition to standard nuclear 
burning. We refer to stars where the energy produced 
by WIMP annihilations is greater than or comparable to 
that from nuclear burning as 'WIMP burners'. Recent 
work has focused upon white dwarfs and neutron stars 
[Ej &\ ■ Here we instead focus on main sequence WIMP 
burners. Previous work in this direction was carried out 
by Salati & Silk some years ago, in the context of 
cosmion dark matter. These authors used an n = 3 poly- 
tropic approximation to estimate the infiuence of WIMP 
annihilation on the main sequence. We instead employ 



full numerical solutions to the hydrostatic equations of 
stellar structure and present results from a simple code, 
as well as supporting preliminary results from a more 
advanced code. We also update the discussion in the 
context of 'modern WIMPs' rather than cosmions, as the 
nuclear scattering cross-sections and the WIMP mass are 
more constrained than in the past. 

The second way WIMPs can influence stellar struc- 
ture is by providing an additional mechanism of heat 
transport in the core. This could reduce the local tem- 
perature gradient [1], potentially inhibiting convection 
and enhancing the pulsation of horizontal branch stars 
d At least for the Sun however, the values cur- 

rentlyfavoured for the WIMP self-annihilation cross sec- 
tion and the upper limit on the WIMP-nucleon cross 
section [2] indicate that this effect is not significant 
In this work, we neglect heat transport by WIMPs, but 
will in a later work include this effect in more detailed, 
time-dependent simulations [T^ . For the current work, 
we have checked that this omission does not affect our 
results significantly. 

Experiments give relatively weak constraints on the 
spin-dependent WIMP-nucleon cross section, so it may 
contribute far more to WIMP capture than the spin- 
independent one. For simplicity, we assume that the 
spin-dependent cross section dominates and is ctsd = 
10~'^*cm^. For a Sun-like star with mass IMq, a cir- 
cular velocity of 220kms~^, a WIMP velocity dispersion 
oiv = 270kms-i and a WIMP mass oflOOGeV, the cap- 
ture rate (as calculated with the full capture expressions 
3 in DarkSUSY ^) is then 



To = 2.90 • lO^^^ s- 



(1) 



We assume a hydrogen mass fraction of 75%, uniformly 
distributed in the star. This is not a good approximation 
for the Sun, but reasonable for newly born stars. The 
capture rate on an arbitrary star is then approximately 
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where is the star's mass and Wcsc its surface escape 
velocity, pw is the ambient WIMP density and TOw is the 
WIMP mass. {(f>) = {v'^sd''^) / ^csc) is ttie average poten- 
tial at the height of the scattering nuclei (hydrogen in this 
case, for which {4>)q — 3.16). For simplicity, we assume 
((/))/ (0)0 = 1. Even if this approximation is not perfect, 
Eq. ([2]) is only necessary for converting between the am- 
bient WIMP density and the capture rate Fc, so it 
is easy to rescale if the reader wishes. We also assume 
V = 270kms~^, though for stars in very different envi- 
ronments to our own this must be adjusted. Throughout 
this paper, we use a WIMP mass of 100 GeV although 
to a first approximation the energy accreted by a star 
immersed in some fixed density of dark matter is only a 
fimction of the cross section, not the mass of the WIMP. 

Once dark matter is captured by the star, it will form 
an approximately thermal (Gaussian) internal distribu- 
tion, of characteristic radius [9,, JJ,, JL4,] 

[ 3fcre 

''w = T^-p; , (3) 

where pc and are the central density and temperature. 
Having been concentrated in the centre of the star, the 
dark matter annihilates with itself, so the equation for 
the evolution of the number of WIMPs in the star over 
time is 



When F(, = 2Fa, the capture and annihilation rates are 
in equilibrium and any dark matter which accretes onto 
the star is instantly converted into additional luminos- 
ity. Here we assume that all of the annihilation products 
interact either electromagnetically or strongly, so that 
they have short mean free paths in the star and the en- 
ergy thermalises quickly. In reality, some fraction of the 
energy will be lost to neutrinos, but this is a small effect 
(typically on the order of 10% or less). The timescale for 
the steady state to be reached will be much less than the 
typical evolutionary timescale of a main sequence star, 
so we assume that this equilibrium has been achieved. 

We studied the effect of this energy input upon stars 
by ahering the FORTRAN code ZAMS [15], which looks 
for stable stellar solutions assuming a constant chemical 
composition and a simple equation of state. The code was 
modified by adding the energy generation due to WIMP 
annihilation to the nuclear energy generation rate. 

Self gravitating systems have negative specific heat, so 
as dark matter injects energy into the core of the star, 
the temperature goes down. Since the nuclear reaction 
rates depend exponentially on temperature, they are sig- 
nificantly reduced for even a small drop in temperature. 
As the ambient density of dark matter is increased and 
the capture rate goes up, the injection of energy due to 
WIMP annihilations eventually reduces the temperature 



enough to shut off nuclear burning. WIMP annihilation 
is then the primary source of the star's luminosity. 

Since WIMP annihilation occurs in a more centralised 
region than nuclear burning, the temperature gradient 
is much steeper in the core of the star than it would 
otherwise be, and the core becomes convective. Outside 
the core, less energy is generated per unit volume than 
if nuclear burning were proceeding normally, so the tem- 
perature gradient is less. The actual temperature is lower 
than in a normal star, but for small dark matter accre- 
tion rates it remains high enough to prevent any major 
increase in opacity, ensuring that energy transport in the 
region above the core remains radiative. The energy from 
the core is easily transmitted through this radiative zone 
to the surface layers of the star. The overall radius of the 
star remains approximately constant while the tempera- 
ture decreases, so the luminosity also decreases. 

As the WIMP accretion rate is raised, the star con- 
tinues to cool until H~ ions are able to survive at in- 
creasingly large depths below the surface. If the star did 
not already have an appreciable surface convection zone, 
one develops. At high enough WIMP capture rates, the 
surface convection zone merges with the inner zone and 
the star becomes completely convective. The addition 
of more WIMPs and hence central luminosity beyond 
this point requires the star to grow in order to transport 
the additional energy to the surface. The two distinct 
stages of this 'evolution' can clearly be seen in Figure 
[TJ This figure shows a black solid line corresponding 
to the zero age main sequence, along with 'evolution- 
ary tracks' plotted as the ambient density of WIMPs is 
increased. These tracks are strongly reminiscent of the 
Hayashi tracks which proto-stars travel along on the fi- 
nal stages of their evolution towards the main sequence. 
The difference is that such stars are shrinking in size 
on the Kelvin-Helmholtz timescale as they radiate away 
gravitational potential energy. In WIMP burners, there 
is constant energy generation in the core (provided the 
WIMP capture rate remains constant) and the stars can 
in principle remain at that position in the HR diagram 
for an arbitrarily long time. 

The analysis of Salati & Silk Q led to qualitatively 
similar conclusions about the backing up of the main se- 
quence along the Hayashi track with increasing WIMP 
density. The shapes of the tracks are different in our 
more detailed analysis, and thus the actual positions of 
the resulting WIMP burners in the HR diagram differ. 
The effects require higher ambient WIMP densities also, 
as we have considered modern WIMPs rather than cos- 
mions (i.e. with lower nuclear scattering cross-sections). 

In Figure [U tracks for stars of mass 0.8 Mq are not 
complete, and show increasingly large gaps for larger 
masses. Using the modified ZAMS code, we could not 
find solutions corresponding to all temperatures between 
the main sequence stars and the cooler Hayashi-like 
WIMP burners. In Figure [2] we plot temperature as a 
function of ambient WIMP density. The discontinuity 
occurs where the temperature of the star drops rapidly. 
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FIG. 2: The temperature of main sequence stars of mass 1 Mq 
( upper curve ) and 0. 7 Mq as a function of WIMP density. A 
spin dependent WIMP-nucleon cross section of a = 10~^* cm^ 
is assumed. 
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FIG. 1: The zero age main sequence of WIMP burners. The 
black solid line is the normal zero age main sequence for stars 
with solar metallicity. The points moving off the main se- 
quence correspond to solutions where stars contain increasing 
amounts of dark matter in their cores. The gaps in the lines 
are addressed in the text. 



shortly before the star becomes fully convective. 

It was not initially clear whether this lack of solutions 
was a real effect or an artefact of the code. We have 
put some effort into understanding whether some of the 
physical simplifications in the code (e.g. surface bound- 
ary conditions, a simplified equation of state and thus too 
rigid a criterion for the onset of surface convection, or the 
treatment of convective heat transport itself) are respon- 
sible for the lack of solutions. We could find no evidence 
to support this idea. Nonetheless, the gaps do not appear 
to have a physical basis. They seem to be a numerical 
artefact arising from the 'shooting' technique used by the 
ZAMS code to solve the boundary-value problem of stel- 
lar structure. In the region of the gaps, small changes 
in the stellar radius cause large changes in the internal 
temperature. The surface defined by the discrepancy be- 
tween the inwardly and outwardly integrated partial solu- 
tions employed in this technique then becomes a highly 
non-trivial function of the model parameters. Finding 
the global minimum of this function is then nearly im- 
possible without an initial guess extremely close to the 
true solution, and the algorithm fails to converge. 

To check that solutions do exist in this region, we inves- 
tigated the gaps using a preliminary version of our next- 
generation WIMP burner code DarkStars. This code has 
been created from generalised versions of capture rou- 



tines in DarkSUSY [iSj which use the full capture rate 
expressions in and the stellar evolution package EZ 
[la l derived from Eggleton's STARS code . The code is 
time-dependent, uses relaxation rather than shooting and 
includes WIMP energy transport, more detailed treat- 
ments of the WIMP distribution, capture rates, equa- 
tion of state, nuclear reaction rates and opacities. Dark- 
Stars and results obtained with it will be described in 
full in an upcoming publication The results of the 
more detailed code (Figure [3]), whilst consistent with the 
ZAMS code close to the main sequence and in the fully 
convective regions on the right-hand side of the HR dia- 
gram, show that intermediate solutions do exist and that 
the gaps are not physical. This cross-check also shows 
that we can trust the solutions of the ZAMS code, but 
shouldn't take the actual gaps in the curves seriously. 

Figure [5] shows that for a spin-dependent WIMP- 
nucleon cross section of cr = lO"-^^ cm^, stars only start to 
change their behaviour when immersed in a dark matter 
density of around 10^ or lO^GeVcm"'^. This is much 
larger than the O.SGeVcm"'^ which is thought to be 
the approximate density of dark matter in the solar sys- 
tem. Most modern simulations of galactic dark matter 
halos do suggest that the density of dark matter should 
be much higher in the centre of the galaxy. This den- 
sity should be made more pronounced by the phenomena 
of adiabatic contraction and dark matter spikes created 
around black holes, though tempered by self-annihilation 
of dark matter and gravitational heating due to the mo- 
tion of stars. It is suggested 18] that densities larger than 
lO^GeVcm"'^ can be found at radii closer than 10"^ pc 
from the galactic centre. Observational challenges asso- 
ciated with constructing an HR diagram from stars in 
this crowded field are severe though, as the very centre 
of the galaxy is shrouded in dust. 
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FIG. 3: Diagram for M« = 2Mq showing that solutions can 
be found corresponding to all temperatures when the more so- 
phisticated DarkStars code is used. 



Finally, we comment upon the 'paradox of youth' im- 
plied by what appear to be hot, young stars close to the 
central black holes of M31 and our own galaxy [l^,[23]- In 
our galaxy, these stars look like ^15 Mq main sequence 
stars with luminosities of the order 1,000 Lq. Such a mas- 
sive main sequence star should only be around 10 Myr 
old, so these objects would appear to have formed very 
recently. However, tidal forces close to the central black 



hole are thought to be too large to be compatible with 
the collapse of a gas cloud to begin star formation. It is 
not easy to make normal UBV measurements of stellar 
temperatures at the galactic centre due to extinction, but 
spectral lines found in the atmosphere of such stars ap- 
pear to be consistent with a temperature of 30,000 K [19l |. 
If this is the case, then it would be difficult to imagine 
that these are normal stars which have been converted 
into fully convective and highly luminous WIMP burners 
from the accretion of dark matter, as such stars would 
have a temperature of only a few thousand degrees. 

A possible explanation for the 'paradox of youth' is a 
combination of collisional stripping of red giant envelopes 
and merger events [l^] , owing to the high stellar densities 
near the centres of the two galaxies. If this is the case, 
then the anomalously hot stars in fact derive from an 
older population. Such a population should therefore also 
contain less luminous, unperturbed stars, including lower 
mass stars still on the main sequence. Provided that 
this explanation is correct, and if main sequence WIMP 
burners exist at all, then such lower mass stars should 
be examples of them. With upcoming observations of 
the galactic centre expected to probe objects as faint as 
1 Lq , a very real possibility exists for the detection of 
WIMP burners in the near future. 
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